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Abstract 

Linear scalar differential equations with distributed delays appear in the study of the local 
stability of nonlinear differential equations with feedback, which are common in biology and 
physics. Negative feedback loops tend to promote oscillation around steady states, and their 
stability depends on the particular shape of the delay distribution. Since in applications the 
mean delay is often the only reliable information available about the distribution, it is desirable 
to find conditions for stability that are independent from the shape of the distribution. We 
show here that the linear equation with distributed delays is asymptotically stable if the 
associated differential equation with a discrete delay of the same mean is asymptotically stable. 
Therefore, distributed delays stabilize negative feedback loops. 



1 Introduction 

The delayed feedback system of the form 

x = F(x,[ [d V (r)}-g(x(t~T),r)), (1) 



is a model paradigm in biology and physics [J, [H, HH, 17, [3, The first argument is the 

instantaneous part and the second one, the delayed or retarded part, which forms a feedback loop. 
The function r] is a cumulative distribution of delays and F and g are nonlinear functions satisfying 
F(0,0) = and g(0,r) = 0. When F : R d x R dxd -> R d and g : R d x R -> R dxd are smooth 
functions, the stability of x — is given by the linearized form, 

y* oo 

x=-Ax- [B(r) ■ drj(T)]x(t -r). (2) 
Jo 

The coefficients A and B(t) <E R dxd are the Jacobian matrices of the instantaneous and the delayed 
parts, r\ : [0, oo) — > R dxd [ s the distribution of delays and (•) is the pointwise matrix multiplication. 
In biological applications, discrete delays in the feedback loop are often used to account for the 
finite time required to perform essential steps before x(t) is affected. This includes maturation and 
growth times needed to reach reproductive age in a population 1131. | la| . signal propagation along 
neuronal axons Q , and post-translational protein modifications I7I. Ildl . Introduction of a discrete 
delay can generate complex dynamics, from limit cycles to chaos [21| . Linear stability properties 
of scalar delayed equations are fairly well characterized. However, lumping intermediate steps into 
a delayed term can produce broad and atypical delay distributions, and it is not clear how that 
affects the stability compared to a discrete delay (lfj| . 

Here, we study the stability of the zero solution of a scalar (d = 1) differential equation with 
distributed delays, 

POO 

x = —ax — b / x(t — T)d/r}(j). (3) 
Jo 
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The solution x(t) £ R is the deviation from the zero steady state of equation (Q]). Coefficients 
a = — Di.F(0, 0) € R and b — — D2-F(0,0) ^ 0, and the integral is taken in the Riemann-Stieltjes 
sense. We assume that 77 is a cumulative probability distribution function, i.e. 77 : R — > [0, 1] is 
nondecreasing, piecewise continuous to the left, tj(t) = for r < and ?7(+oo) = 1. Additionally, 
we assume that there exists v > such that 



e vr dt]{r) < 00. (4) 



This last condition implies that the mean delay value is finite, 

/>OC 

E = I Tdr}(r) < 00. 
Jo 

The corresponding probability density function is /(r) given by dr)(r) — f(r)dr, where the deriva- 
tive is taken in the generalized sense. The distribution can be continuous, discrete, or a mixture of 
continuous and discrete elements. When it is a single discrete delay (a Dirac mass), the asymptotic 
stability of the zero solution of equation ^ is fully determined by the following theorem, due to 
Hayes [l2], 

Theorem 1. Let /(r) = 5(t — E) a Dirac mass at E. The trivial solution of equation {SJj is 
asymptotically stable if and only if a > \b\, or if b > \a\ and 

E < arccos(-a/b) ^ 



y/b 2 - a? 

There is a Hopf point if the characteristic equation of equation (|3|) has a pair of imaginary 
roots and all other roots are negative. For a discrete delay, the Hopf point occurs when equality 
in (f5|) is satisfied. Moreover, for any distribution 77, there is a zero root along the line —a = b. At 
—a = b = 1/E, there is a double zero root. When a > —1/E, all other roots have negative real 
parts, but when a < —1/E, there is one positive real root. Thus, the stability depends on 77 if and 
only if b > \a\. Moreover, only a Hopf point can occur when b > |o|. Therefore, a distribution 
of delays can only destabilize equation (|3|) through a Hopf point, and only when b > \a\. This is 
a common situation when the feedback acts negatively on the system (D 2 F(0, 0) < 0) to cause 
oscillations. 

Assuming b > and making the change of timescale t — > bt, we have a — * a/6, b — > 1 and 
7](t) — > T)(br). Equation (J3J) can be rewritten as 

o 

x(t- T)dri(T). (6) 







The aim of this paper is to study the effect of delay distributions on the stability of the trivial 
solution of equation ©, therefore, we focus on the region \a\ < 1. To emphasize the relation 
between the stability and the delay distribution, we will say that n (or /) is stable if the trivial 
solution of equation © is stable, and that 77 (or /) is unstable if the trivial solution is unstable. 

It has been conjectured that among distributions with a given mean E, the discrete delay is 
the least stable one [H, @. If this were true, according to Theorem [TJ all distributions would be 
stable provided that 

„ arccos(— a) , . 

E < 7=^r - (7) 



This conjecture has been proved for a = using Lyapunov-Razumikhin functio ns 1X511 , and for 
distributions that are symmetric about their means [f(E — r) = f(E + r)] [4|, Is LLJ, UM- ^ has 
been observed that in general, a greater relative variance provides a greater stability, a property 
linked to geometrical features of the delay distribution 0. There are, however, counter-examples 
to this principle, and there is no proof that for a 7^ the least stable distribution is the single 
discrete delay. It is possible to lump the non-delayed term into the delay distribution using the 
condition found in [15j . but the resulting stability condition, Ej(\ + a) < 7r/2, is not optimal. 



2 



Here, we show that if inequality ([7]) holds, all distributions are asymptotically stable. That is, 
distributed delays stabilize negative feedback loops. 

In section [21 we set the stage for the main stability results. In section [31 we show the stability 
for distributions of discrete delays. In section [H we present the generalization to any distributions 
and in section 03 we provide illustrative examples. 



2 Preliminary results 

Let T) be a distribution with mean 1. We consider the family of distributions 

frjir/E), E>0, 
VEiT) = \H(T), E = 0. (8) 

where H(t) is the step or heaviside function at 0. The distribution tje has a mean E > 0. The 
characteristic equation of equation ([6|), obtained by making the ansatz x(t) = cxp(— Xt), is 

poo 

X + a + e- XT dr] E (T) = 0. (9) 
Jo 

When condition ^ is satisfied, the distribution ije is asymptotically stable if and only if all roots 



of the characteristic equation have a negative real part Re(A) < [22} . Condition guarantees 
that there is no sequence of roots with real parts converging to a non-negative value. The leading 
roots of the characterisitc equations are therefore well defined. When E = 0, i.e. when there is no 
delay, there is only one root, A < 0. When E > 0, the characteristic equation has pure imaginary 
roots A = ±iu only if < lu < lu c — i/l — a?. Thus, the search for the boundary of stability can 
be restricted to imaginary parts lu € (0,u> c ] [6|. 
We define 

/•OO 

C(w) = / cos(u;t)*7e(t), (10) 







oc 



S(u) = / sin(wr)d77 B (T). (11) 
Jo 

We use a geometric argument to bound the roots of the characteric equation of equation © by 
the roots of the characteristic equation with a discrete delay. More precisely, if the leading roots 
associated to the discrete delay are a pair of imaginary roots, then all the roots associated to the 
distribution of delays have negative real parts. We first state a criterion for stability: if S(lu) < ui 
whenever C(lu) + a = 0, then / is stable. The larger the value of S(u>), the more "unstable" the 
distribution is. We then show that a distribution of n discrete delays /„ is more stable than an 
certain distribution with two delays /*, i.e. S n (uj) < S*(lo). We construct /* and determine that 
one of the delays of this "most unstable" distribution /* is rj* = 0, making it easy to determine 
its stability using Theorem [Q We then generalize for any distribution of delays. 

The next proposition provides a necessary condition for instability. It is a direct consequence 
of theorem 2.19 in [22] . We give a short proof for completeness. 

Proposition 1. If the distribution r\E is asymptotically unstable, then there exists uj s G [0, £J C ] 
such that C{u) s ) + a = and S{lo s ) > lo s . 

Proof. Suppose that the distribution t\e is asymptotically unstable. The roots of the characteristic 
equation depend continuously on the parameter E and cannot appear in the right half complex 
plane. Thus there is a critical value < p < 1 at which r/ p E loses its stability, and this happens when 
the characteristic equation (JH]) has a pair of imaginary roots A = ±iw with < lu < uj c = yl — a 2 , 
i.e. through a Hopf point. Splitting the characteristic equation in real and imaginary parts, we 
have 

/•OO 

cos(ujT)dr]pE (t) + a = 0, 
sm(u)T)dr) P E(T) = lu. 
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Rewriting in term of tje, we obtain 

"OC 

cos(ujpT)dr]E(T) + a = 0, 

sm(u!pT)drjE(r) = lo. 

Finally, denoting lo s = pto < lo < lo c , we have 

cos(u> s T)dr]E(T) + a = 0, 

sin(o; s T)dr] E {T) = uj = uj s /p > uj s . 

This completes the proof. □ 

Proposition [T] provides a sufficient condition for stability: 

Corollary 1. The distribution rjE is asymptotically stable if (i) C(uj) > —a for all lo £ [0,o; c ] or 
if (ii) C(lo) = —a, lo G [0,o; c ], implies that S(u>) < lo. 

Proposition [T] suggests that the scaling r/E — t](t/E) is appropriate for looking at the stability 
with respect to the mean delay. The mean delay scales linearly, and unstable distributions therefore 
lose their stability at a smaller values of the mean delay, under this scaling. The condition S(lo s ) < 
lo s is however not necesssary for stability, as one can find cases where S{lo s ) > lo s even though the 
distribution is stable. This happens when an unstable distribution switches back to stability as E 
is further increased (see for instance or [j| and example I5.3p . 

3 Stability of a distribution of discrete delays 

We define a density of n discrete delays Tj > 0, and pi > 0, i = 1, n, n> 1, as 

71 

Mr) =J2pAr-Ti) (12) 

z=l 

where S(t — Ti) is a Dirac mass at Tj, and 

n n 

*^2piTi = E, and ^^Pi = 1. 

In this section, we show that /„ is more stable than a single discrete delay. We do that by observing 
that among all n-delay distributions, n > 2, that satisfy C n (uj s ) + a = for a fixed value lo s < oj c , 
the distribution /* that maximizes S n (w s ), 

m&x{S n (Lo s )\C n (Lo s ) + a = 0} = S*(lo s ), (13) 

has 2 delays. We show that S*(lo s ) < lo s , implying that all distributions are stable. The following 
lemma shows how to maximize S(lo s ) for distributions of two delays. 

Lemma 1. Let fi be a delay density with mean E. Assume in addition that there exists lo s € (0, lo c ) 
such that C(to s ) = —a < cos(lo c E). Then there exists t£, pi and P2 such that 

rt - 0, (14) 
Pl+P* 2 = l, (15) 

p\ +p* 2 COs(w s T 2 *) = pi COs(lO s Ti) +p 2 COs(w s T 2 ), (16) 

p* 2 r* = E. (17) 



4 



Moreover, there is at most two solutions for t% with r| < tt/uj s . If t£ is the smallest solution, we 
have that r| < t 2 and 

2 

S*(u> s ) = sin ( w s T *) > S{lo s ). 

i=l 

Proof. To see that there is always a solution, let c > be the smallest value such that the inequality 
cos(#) > 1 — c9 is verified for all 9. [c = 0.725... by solving c = sin(0), with 1 — 9sm(9) — cos(#).] 
We have that 1 — cojE < C(lo) < cos(u>E). Thus, the line 1 — dcuE that goes through C(uj s ) satisfies 
d = (1— C(u} s ))/ (w s E) < c, and therefore crosses the curve cos(ui s E) at some points. The smallest 
solution t"! is the one such that 1 — dw s r| = cos(w s T2 ). This way, 

Pl +p* 2 cos(cj s t*) = 1 - (1 - C(u.))pZt2/E, 
= C(lo s ). 

These new delay values maximize S(u s ) under the constraints that C{u> s ) + a = and that the 
mean remains E. That is, we will prove that 

2 2 
S*(w„) = 22 Pi siri ( w sT*) >^Pi sin(w s Ti), 

»=1 i=l 

for all admissible pi, Ti. Two show that, we recast the problem in a slightly different way. Writing 
u = lo s t\, v = to s T2 and T = lo s E, we can express parameters pi in terms of (u, v): 

v — T T — u 

Pi = and p 2 = , 

V — u v — u 

where u < T < v. We consider C and S as functions of (u, v). 

v — T T — u 

C(u, v) = cos(w) H cos(u), (18) 

v — u v — u 

v — T T — u 

S(u, v) = sinful H sin(w). (19) 

v — u v — u 

Equation (|19p is to be maximized for (u, v) along the curve h = {u, v} implicitly defined by 
the level curves C(u,v) = —a. There are either two solutions in v, including multiplicity, of 
the equation C(0,v) — —a or none, so the curve can be parametrized in a way that («(£),«(£)) 
satisfies (it(0) = 0, v(0) = v max ) and (u(l) = 0, v(l) = v min ), with v min < v max . We claim that S 
is maximized for £ = 1, i.e. u — and v — v max . This is true only if S{u{^), v(£)) is increasing 
with ^. That is, the curve h must cross the level curves of S upward. It is clear that S is a 
decreasing function of v, for u fixed and an increasing function of u, for v fixed. Thus the level 
curves S(u,v) — k can be expressed as an increasing function vs,k(u) such that 

S(u,vs,k(u)) = k, 

when k is in the image of S. Likewise, equation Q18[) can be solved locally to yield vc, a (u) such 
that 

C(u,v c ,a(u)) = -a-, 

whenever —a is in the image of C. The function vc, a {u) could take two values on the domain of 
definition. Because S is decreasing in v, we choose the lower solution branch for vc, a (u)- If, along 
that lower branch, the slope of vc, a {u) is larger than that of vs,k{u), then as v decreases along the 
curve c, S increases. Therefore, to show that (0,t>o(0) = u m in) maximizes S, we need to show 
that 

dvcA u ) > dvs,k( u ) > q ^ 20 ^ 
du du 
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Figure 1: How delays are replaced to get an maximal value of S* . In this example, a = —0.1. Parameters are 
u = 0.2, v = 2, pi = 0.37, p 2 = 0.63 and T = 1.33. The parameters maximizing S are u* = 0, v* = 1.76, pi = 0.24 
and P2 = 0.76. 



It is clear that dv s ,k{u) / du > 0. The pointwise derivatives of the level curves at (u,v) are 

dvc{u) v — T — cos(u) + cos(u) + (v — u) sin(w) 

du T — u cos(u) — cos(u) — (v — u) sin(w) ' 

dvs{u) v — T — sin(u) + sin(u) — (v — u) cos(u) 

du T — u sin(u) — sin(w) + (v — u) cos(u) 

Because only the lower branch of vc is considered, we restrict (u, v) where dvciu) / du < +oo. This 
is done without loss of generality since S is strictly larger on the lower branch than on the upper 
branch. Along the lower branch, vc(u) < it. Inequality PU|) then holds if 

(v — u)[2 — 2(cos(m) cos(v) + sin(w) sin(v)) + (v — u)(sin(u) cos(v) — cos(u) sin(u))] > 0. 

Notice that this inequality does not depend on T, which cancels out, nor on a, since comparison 
is made pointwise, for any level curves. The inequality can be simplified and rewritten in terms of 

z = v — u > 0, 

z [2 - 2 cos(z) - z sin(z)] > 0. 

It can be verified that this inequality is satisfied for z e (0,7r]. Therefore, S is maximized when 
u = and v = v c , a (0) = v min . □ 



Theorem 2. Let f n be a density with n > 1 discrete delays and mean E satisfying inequality 
The density f n is asymptotically stable. 

Proof. Single delay distributions (n — 1) are asymptotically stable by Theorem [T] We first show 
the case n = 2. 

Consider a density fa, with t± < t%. Suppose C{oj s ) + a = for a value of uj s < lo c (if 
not, Corollary Q] states that fa is stable). Remark that —a = C(uj s ) < cos(uj s E). Indeed, from 
inequality J7J) and to s < u> c = \[\ — a 2 , we have cos(oj s E) > cos(uj c E) > —a. Replace the two 



G 



delays by two new delays with new weights: rj* = and t£ > the smallest delay such that the 
following equations are satisfied: 

P*2 T 2 = Pi T i +P2T2, (21) 
pi + p* 2 cos(w s T2 ) = pi cos(w s ti) + P2 cos(w s r 2 ), (22) 
P*l+P*2=Pl+P2 (-1). (23) 

Lemma [T] ensures that there always exists a solution when C(oj s ) < cos(lo s E). Additionally, r| < r 2 
and 



S*(uj s ) = },p* sin(w s r*) > ^Pi sm(ui s T. t ). 



i=i i=i 

That is, the new distribution * maximizes the value of S. Therefore, if we are able to show 
that distributions with a zero and a nonzero delay satisfy S(u> s ) < ui s , then by Corollary [I] all 
distributions with two delays are stable. Consider /(r) = (1 — p)S(t) + pS(r — r). Suppose that 
there is u> s < ui c such that 

C(lu s ) = 1 — p + pcos(uj s r) = —a. 

We must show that S(lu s ) = psin(o; s ) < lu s . Summing up the squares of the cosine and the sine, 
we obtain 

p 2 = {-a+p-\f + S 2 (u Js ), 

so 



= vV-(-a+p-l) 2 - 
By assumption, the mean delay statisfies inequality (J7J, 

arccos(— a) 
pr < — - 

Thus, 

arccos(— (a + 1 — p)p~ 1 ) , arccos(— (a + 1 — p)p~ 1 ) 

uj s = > pVl - a ; ; • 

r arccos(— a) 

Because (a + 1 — p)/p > a for p 6 (0, 1] and a 6 (—1, 1), we have the following inequality 

arccos(— a) ^ arccos(— (a + 1 — p)p^ 1 ) 



' ^l-((a + l-p)p^) 2 

Thus, 

at ^ n 1 T T\2 s R 2 - arccos(-(a + l-p)p- 1 ) 

S(u) s ) = VP - {-a+p- I) 2 < pVl - a 2 * — r - < lu s . 

arccos(— a) 

This completes the proof for the case n = 2. 

For distributions /„ with n > 2 delays, the strategy is also to find a stable distribution that 
keeps C{(jJ s ) constant and increases S(uj s ), assuming that C(uj s ) + a = 0. This requires two steps. 
In the first one, all pairs of delays Ti < Tj for which the inequality 

^ PfcCOs(w s T fc ) < COsf UJ S ^ PkTk), (24) 

ke{i,j} ^ ke{i,j} ' 

holds are iteratively replaced by new delays r* = and r* < Tj , as done in Lemma [T] This 
transformation preserves E, C(uj s ) and increases S(iv s ). This is repeated until there remains 
m < n delays with n > 0, i = 2, m such that 



^ Pk cos(uj s T k ) > cosf lo s ^2 PkT k Y 

ke{i,j} ^ ke{i,j] ' 
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Figure 2: Stability chart of distributions of delay in the (a/b, bE) plane. The distribution-independent stability 
region is to the right of the blue curve. The distribution-dependent stability region is the shaded area. All stability 
curves leave from the point (a = —b, E = 1/6). The signs of the real roots of the characteristic equation Aq, Ai along 
a = — b are distribution-independent. 



for i/je {2, m}, and t\ — 0. (The rj are not the same as in the original distribution, the * 
have been dropped for ease of reading.) The positive delays Tj > satisfy 

m m 
y^ y PjCOs(uJ s Ti) > COS (m s y^PiTj). 

i=2 i=2 

while, by assumption, 

m 

^^Pi cos(u> s Ti) = —a < cos(u> s E). 

i=i 

The second step is to replace all delays r i; i = 2, ...,m with the single delay f 2 = Y^i=2"Pi T i- 
We now have a 2-delay distribution with f i = and T2 > 0, pifi +P2J2 = E, C(uj s ) < C(uj s ) and 
S(u) s ) > S(lu s ). Replace T2 by the delay < r 2 so that C*(uj s ) — —a, while keeping E constant. 
Existence of is shown using the notation from the proof of Lemma [lj and noting that C(0, v) 
and S(0, v) are both decreasing in v. This change of delay has the effect of increasing S: S*(uj s ) > 
S(uj s ). Therefore, we have found a pair of discrete delays (0,r|) such that C*(lu s ) = C(lu s ) and 
uj s > S*(u! s ) > S(lu s ). By Corollary [TJ /„ is asymptotically stable. □ 



4 Stability of a general distribution of delays 

From the stability of distributions of discrete delays to the stability of general distributions of 
delays, there is a small step. First we need to bound the roots of the characteristic equation for 
general distributed delays. 

Lemma 2. Let T]e be a delay distribution with mean E satisfying inequality There exists 

a sequence {rj nt E}n>i with distribution r\ n ,E having n delays, such that rj ni E converges weakly to 
rjE- Then X is a root of the characteristic equation if and only if there exists a sequence of roots 
A n for r\ n ,E such that lim„^oo A„ = A. Let {^ n } n >i be a sequence of real parts of roots of the 
characteristic equations. Additionally, 

limsup/i„ = fj, < 0. 
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Proof. Consider A„ = fi n + ioj n a root the characteristic equation for rj n . E . E satisfies inequality 
0, so fj, n < 0. So 



A„ + a + l t 
o 



x - T d[rf E {T)-rf n>E {r)] 



-A„t 



dr)„ tE (r) 



e- x - T d[i lE (r) -T) n , E (r)] 



0, 



as n — > oo by weak convergence. Thus any converging sub-sequence of roots converges to a root 
for rj E . The same way, if A is a root for r] E , 



X + a+ e~ AT dr] niE {T) 



A + a + / e- XT d[ri^ E {T) ~ t] E {t)} + / e~ XT d m { T ) 
o Jo 



-Ar 



d[»?n,E(r) - J?s(r)] 



0, 



as n — > oo. Convergence is guarantedd by inequality ([!]). Therefore, each root A lies close to a 
corresponding root A„. 

Denote fi — \iu\s\xp n _ too fi n . Then /i is the real part of a root of the characteristic equation 
associated to r/ E . fi n < for all n, so /i < 0. Suppose fi = 0. Without loss of generality, we can 
assume that all other roots have negative real parts. Then r\ E is at a Hopf point, i.e. the leading 
roots of the charateristic equation are pure imaginary. Consider the distribution r\a.p(j) — v( T /p) 
and the associated real parts fia.p, where the subscript a is there to emphasize the dependence of 
the stability on a. Then, by continuity, there exists (a,p) in the neighborhood e > of (a,E) for 
which r/a.p is unstable, i.e. fia,p > 0. For sufficiently small e > 0, inequality ((7j) is still satisfied: 



arccos(— a) 



Additionally, r\ n ^ p converges weakly to rj a . p - However, because ?/ a ,p is unstable, there exists N > 1 
such that r] n ,a,p is unstable for all n > N, a contradiction to Theorem [5] Therefore fi < 0. □ 

Theorem 3. Let rj E be a delay distribution with mean E satisfying inequality UN, The distribution 
rj E is asymptotically stable. 

Proof. Consider the sequence of distributions with n delays {ri nyE } n >i where rjn iE converges weakly 
to r\ E . By Lemma the leading roots of the characteristic equation of r\ E have negative real parts. 
Therefore r/ E is asymptotically stable. □ 

Is there a result similar to Theorem |31 for the most stable distribution? That is, is there a 
mean delay value such that all distributions having a larger mean are unstable? When a > 0, 
the answer is no. For instance the exponential distribution with parameter is asymptotically 
stable for all mean delays, a property called unconditional stability. Other distributions are also 
unconditionally stable for a > 0. Anderson has shown that all distributions with smooth enough 
convex density functions are unconditionally stable 0, but densities do not need to be convex to 
be unconditionally stable. For example, the non-convex density /(r) = 0.5[<5(r) + 8(t — 2E)] has 
mean E but is unconditionally stable. However, no distribution is unconditionally stable for all 
values of a £ [—1,0), although some are for a > a* with a* > — 1 (see example below). 

From the results obtained here, we have the most complete picture of the stability of equation 
when the only information about the distribution of delays is the mean (figure \2§ ■ 
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Corollary 2. The zero solution of equation {3|) is asymptotically stable if a > —6 and a > |6| or 
ifb> \a\ and 

„ arccosf— alb) 

E < ; - , 

Vb 2 - a 2 

The zero solution of equation^ may be asymptotically stable (depending on the particular distribu- 
tion) if b > \a\ and 

arccosf— alb) 

E > , . 

~ Vb 2 - a 2 

The zero solution of equation (0) is unstable if a < —b. 

5 Boundary of stability 

The exact boundary of the stability region in the (a, E) plane can be calculated by parametrizing 
(a(u), E(u)). Consider the distribution n. Then, at the boundary of stability 

oo 



= iuj + a+ / e~ tu > T dr](T/E), 
Jo 

/•OO 

= ioj + a+ / e~ lu)ET dr]{T), 



setting u — Et, 



" 



poo 

= i- + a+ / e~ iUT drj{T) 







E 

Separating the imaginary and the real part, we obtain 

ii 

a(u) = -C(u) and E(u) = (25) 

S{u) 

for u > 0. The fact that u depends on E is not a problem: u — > oo if and only if E — > oo, and 
u — > if and only if E — > 0. Equations (|^|) allows systematic exploration of the boundary of 
stability in the (a, E) plane. 

5.1 Exponential distribution 

The exponential distribution /(t) = e~ r has normalized mean 1, and 

1 u 

= T~i 2 aild S ( U ) = 1 j_ 2 " 

1 + 1 + IT 

The stability boundary is given by E = — 1/a, for —1 < a < 0. Therefore the exponential 
distribution is not unconditionally stable for a < 0. 

5.2 Discrete delays 

The exponential distribution is also not the most stable distribution. The density with a zero and 
a positive delay is /(r) = (1 — p)S(r) +pd(r — r), p 6 (0, 1]. After lumping the zero delay into the 
undelayed part, the exact stabiltity boundary becomes 



arccos ( — (a + 1 — p)p x ] 

E = pr = 



'l-((a+l-p)p-iy 
This has an asymptote at a — 2p — 1, which can be located anywhere in (—1, 1] 
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In general, for a distribution with n delays, 

n 

a(u) — — ^y^^Pi cos(nTj) and E(u) 



E"=iP» sin(uTi)' 



The boundary of the stability region can be formed of many branches, as with a distribution with 
three delays in figure [31 

5.3 Gamma distribution 

As the mean E is increased, distributions can revert to stability. This is the case with the second 
order gamma distribution (also called strong kernel) with normalized mean 1, 

f(r) = 2 2 re- 2 \ (26) 

We have 

C(u) = 7T, and S(u) — 

(l + « 2 ) 2 (1 + u 2 ) 2 

The boundary of stability is given by 



(«w^(")) = (^V^' (i+u)2 )- 



There is a largest value a = 0.1216. For large values of E, a — > + . Therefore the boundary of 
the stability region is not monotonous; for a £ (0, a), / first becomes unstable and then reverts to 
stability as the mean is increased (figure [3j. 
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